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ABSTRACT 



Previous analysis about the deparametrization and path integral quantization of cos- 
mological models are extended to models which do not admit an intrinsic time. The 
formal expression for the transition amplitude is written down for the Taub anisotropic 
universe with a clear notion of time. The relation existing between the deparametrization 
associated to gauge fixation required in the path integral approach and the procedure of 
reduction of the Wheeler-De Witt equation is also studied. 
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1 Introduction 



A time dependent canonical transformation can turn a parametrized system into an or- 
dinary gauge system. Indeed, in reference it has been shown how to build a canonical 
transformation such that the fixation of the new coordinates is equivalent to the fixation 
of the original ones. More recently, and within this context, simple cosmo logical models 
were used as examples to show that gravitation can be quantized as an ordinary gauge 
system if the Hamilton- Jacobi is separable. In this paper we discuss how to extend this 
procedure in order to study the case of Hamiltonian constraints which do not allow for 
an intrinsic global time We treat the case of the Taub anisotropic universe as an 
example of physical interest with the intention to present a consistent deparametrization 
procedure which leads to write down a formal expression for the quantum transition am- 
plitude. We study the realization of the analysis of the reduction within the context of 
the path integral approach and we show that the extrinsic time identified in || coincides 
with one which can be obtained by means of the systematic method proposed previously 
in references We will also discuss that this coincidence is, in fact, an evidence 

of the narrow relation existing between the canonical approach of quantization and the 
path integral quantization with a clear notion of time. Then, we analyse aspects of the 
structure of the constraint surface in order to establish the correspondence. 

In section 2 we analyse the steps of deparametrization and path integral quantization 
of the Taub universe. Then, in section 3, we discuss the canonical quantization, leaving 
the conclusions for section 4. 

2 Path integral approach 

Let us consider a Hamiltonian constraint H with the generic form 

N 

H = G^TTiTTj + C ne a?91 ~ (1) 
n=0 

3 



where G 1 ^ represents the minisuperspace reduction of a supermetric with Lorentzian sig- 
nature and {g l , n 1 } are the degrees of freedom of the system and its conjugate momenta 
respectively. The indices i,j G {0, ...d} and n G {0, ...N}, being N < d. 

In fact, it is possible to prove that if the matrix with components {a™} admits a 
dimensional extension {a;™}— >{af} leading to obtain a matrix a\ of dimension d + 1 
that is diagonalizable in a base {x 1 } in which the form of the supermetric Cry is also 
diagonal, then a canonical transformation exists leading to obtain a classically equivalent 
Hamiltonian constraint which admits an intrinsic global phase time. 

An example of this kind of constraint (|l|) is given by the Taub cosmological model, 
which represents a particular case of the anisotropic generalization of the Friedmann- 
Robert son- Walker universe with curvature k = +1 (i.e. by setting to zero one of the 
degrees of freedom of the diagonal Bianchi type- IX universe). Other models of this class 
can be found in reference |4]]. In this paper, we explicitly apply the procedure described 
above to the case of the Taub universe in order to turn this cosmological model into an 
ordinary gauge system and impose on it canonical gauges to identify a global phase time. 

The action functional of the Taub universe reads 

with the Hamiltonian constraint 

H = e- m (ix 2 + - Til) + l -e n {e-^ - 4e~ 2 ^) « 0. (3) 
Because the factor e _3Q is positive definite, this constraint is equivalent to 

H = n l _ ^ + I e 4n (e -8/3 + _ 4e -2/3 +) ^ q_ (4) 

3 

Here, we identify the previous nomenclature with physical parameters as (g°, (yf 1 ) = (f2, (3 + ). 

Besides a separable Hamiltonian, the Taub cosmological model includes true degrees 
of freedom and a potential which vanishes at a given point of the phase space, so making 
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impossible the definition of an intrinsic time in terms of the original variables. As it is 
easy to see, for (3 + = — | In 4 the potential is zero. 

Note that in the case (§]) the parameters of the Hamiltonian are given by a\ = a\ = 4, 
a\ = —8, a 2 = —2 and c\ = — 4c2 = 3. The Hamiltonian is not separable in terms of 
the coordinates and momenta (ft, /3 + , 7Tq, tt + ). Then, in order to apply the method of 
deparametrization we define the coordinates 

x = tt - 2f3 + , y = 2tt-f3 + (5) 

and, thus, we can write 

11 = ^-^ + 1^ -4e 2 ^0. (6) 

At this stage we then have H = Hi(x,it x ) + H 2 (y, 7r y ) with Hi > and H 2 < 0. 
Since the potential vanishes for y = 2x — In 2, a time in terms of only x, y does not exist. 
Hence before turning the model into an ordinary gauge system we shall make a first 
canonical transformation to the variables (x, s,it x ,tt s ) so that the new potential has only 
one positive definite term. The resulting new coordinates will correspond to a set {q 1 } in 
terms of which an intrinsic time exists. 

We shall perform a canonical transformation matching H 2 (y, 7r y ) = — 7rj?, so that ir s = 



±y— H 2 (y, iTy). This is achieved [[| by introducing the generating fuctions of the first 
kind 
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fi(y,s) = ±-e y sinhs. (7) 



The momenta are then given by 



so that 



ix y = i-e^sinhs 
o 

7r s = ±-e y coshs 
3 



S 2 +^=^(sinh 2 S + l) = 7 r s 2 



and the Hamiltonian can be written as 

H(8, X, 7T„ 7T X ) = -7T 2 S + 7T* + ^ « 0. (9) 

It is important to note the reason why we have introduced two possible definitions of 
/i, indeed this ambiguity is reflected by the fact that the constraint can be written as 



H = [ -7T S + y 71-2 + Igtoj |^ + ^2 + l e 4x ] ~ 0) (10) 

and if we had chosen a definite sign, according to the resulting sign of tt s , only one of 
the factors would be zero. But we shall consider the Hamiltonian constraint (|9|) as the 
starting point for applying the deparametrization procedure, and we shall not go back 
to put the results in terms of the original variables. Indeed, at the level of the variables 
(x, s, ii x , ii s ) there is no justification to prefer one possible sign of tc s . 

A point to be noted is the fact that the action S[x, s, ir x , ir s , N] will differ from 
S[fl, P + , tcq, 7t + , N] in surface terms associated to the transformation generated by fi, 
namely D: 

S[x, s, tt x , 7r s , N] = J ^ x ^ + Tr 8 ^-NH(x,8,n x ,n B )jdT (11) 
= S[n,P + ,n a ,n+,N] + [D(T)]Z, (12) 

and when we turn the system into an ordinary gauge one we will have additional 
surface terms B(t); but as we will not be interested in a transition amplitude between 
states labeled by the original coordinates {x, y}, but by {x, s}, we shall not require that 
D(t) + B(t) = 0, but only that B(t) = in a globally good gauge such that an intrinsic 
time is defined JIJ . This ensures that the action of the ordinary gauge system weights the 
paths in the same way that the original parametrized action. 
Now, let us introduce the coordinates 

u = L e z(*+s) v = L e *(*->) (13) 

12 ' 12 K ' 



which lead to the equivalent constraint 

H' = n u 7T v + 1. (14) 

After two succesive canonical transformations in order to obtain a system with non 
conserved observables UH we find the following definition of variables of the gauge 
system, which arises from indentifying P with the energy and going back to (x, s, n x , ir s ): 

p 2(x-s) 

^ 12P ' 

Q = ^e 2 ^) + -l(le 2 (— )(l-P )-r ? T(r)), 

= 9(n 2 x - 7r 2 )e- 4 * + 1, 

P = 3(7r s + 7r x )e- 2 ^ +s l (15) 

Here T(r) is a monotonous function of r and rj = ±1. By fixing the gauge \ = Q°—T(t) = 
it is possible to see that an extrinsic time is 

1 e Ax 

t(x, s, ir x , n s ) = — ■ . (16) 

36 7T X + 7T S 

The constraint surface splits into two sheets given by the sign of ir s . On each sheet 
the intrinsic time can be defined as 

t(x,s) = ^gn(n s )e 2 ^- s \ (17) 

which is associated to the canonical gauge x = Q°P — 7 7^( r ) = because P is propor- 
tional to 7i s + n x and sign(ir s + tc x ) = sign(7T s ). The end point terms associated to the 
transformation (x, s, ir x , ir s ) — > (Q\Pi) are |] 

B{t) = 2Q -Q P -2 V T[T ' ! 



P 

1 f \ e ^ _ 2r ] e 2{x+s) T( 



3(vr s + IX x ) V6 

and with the gauge choice defining an intrinsic time they vanish on the constraint surface 
Pq = 0. The expression for the quantum propagator is 



(x 2 , s 2 \x-i, St) — J DQDPexp i (pdQ - -^dT 



(19) 
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The end points are given as T\ = i^e 2 ^ 1-81 ' and T2 = ±j^e 2 ^ X2 ~ S2 \ Note that in the 
gauge defining the intrinsic time the new coordinate Q coincides with j^e 2<yX+s \ so that 
the paths go from Q 1 = j^e 2 ( Xl+Sl ' to Q 2 = j^e 2< ^ X2+S2 \ We have obtained a propagator 
with a clear distinction between time and the physical degree of freedom. The path 
integral corresponds to that for a conservative system with Hamiltonian rj/P. Because 
1] = sign(7r s ) then at level of the physical degrees of freedom we have two disjoint theories, 
one for each sheet of the constraint surface. 

Propagators of the form ([19]) have been previously treated within the context of the 
minisuperspace approximation to quantization of homogeneous cosmological models. In 
fact, the application of the deparametrization procedure can be extended to several sys- 
tems such as the case of the Kantowski-Sachs model, the homogeneous Szekeres universe, 
dilaton cosmological models and the closed de Sitter model (see [f| and references therein 
for a complete review). 

On the other hand, a global phase time can also be identified among the coordinates 
labeling the quantum states. Because the momentum ir s does not vanish on the constraint 
surface, the coordinate s is itself a time. More precisely, on each sheet of the constraint 
surface we can define an intrinsic time in the form 



because [t*,H] = 2w s sign(ir s ) > 0. Note that this time can be directly recognized from 
eq. ([TOD , as it was shown in reference 0. 

Although we do not use this time as the time parameter in the path integral, the 
interpretation of the result can be made more clear by recalling that one of the coordinates 
which identify the states is a global phase time. In fact, we can write the transition 
amplitude as 



t* = —s sign(n s ) 



(20) 




(21) 
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With this interpretation, the propagator becomes completely analogous to that of a 
mechanical system with a physical degree of freedom x whose evolution is given in terms 
of a true time t*. 

At this point it is natural to ask wether this time could have been obtained in a direct 
way with our deparametrization procedure. This is in fact possible, even in the case that 
we include a matter field in the model. Consider the Hamiltonian constraint for the Taub 
universe with a massless non interacting scalar field (a massless dust): 

H = -7v 2 n + ttJ + n 2 + + ie 4n (e- 8/3 + - Ae~ 2 ^) « 0. (22) 

If we change to the coordinates x and y, and then perform the canonical transformation 
generated by fi(y, s), we obtain the equivalent constraint 

H = -tt 2 + + V » 0, (23) 

where we have redefined tc^ — > tt^/^/S. The corresponding Hamilton-Jacobi equation is 
separable: 

and the solution is clearly of the form W\{x, ir x ) + W 2 (<f>, vr^) + W 3 (s, ir s ). Introducing the 
integration constants b 2 = ix 2 , and a 2 such that a 2 + b 2 — E = n 2 we obtain 



W = sign(ir x ) J dx]ja 2 — -e 4x 

+ s sign(TT s )Va 2 + b 2 — E + sign^^VtP. (25) 



If we match E = P Q then we have 



dP 



Po=0 



-siom-Ks) — . = . 



As [CT,Po\ = 1 then we can inmediately define an extrinsic time as 



t( s ,7r s ) = Q° = ~- (26) 



As in the variables (s, x, 0, 7r s , ir x , n^) the constraint surface is topologically equivalent to 
two disjoint half planes, each one corresponding to 7r s > and to tt s < 0, we can define 
the time as 



which coincides with the time t* found just before following the steps of the reduction 
proposed in reference ||. In terms of the new variables this is an intrinsic time; of course, 
in terms of the original variables it is an extrinsic time. 

All the results include the sign of the momentum tt s , which comes from the double 
sign in the definition of f\; this fact manifestly shows which is the counterpart of the 
ambiguous sign in the definition (|j) within the context of the path integral approach. 
However, some authors have suggested that the constraint of a parametrized system, which 
is linear in the momentum conjugated to the time, may be hidden in the Hamiltonian 
formalism for the gravitational field. According to this point of view, one should choose 
only one of both possible signs for the generator fx , and there would not be two coexisting 
theories (see || for an interesting discussion about this point). Our results, instead, reflect 
that we consider the quadratic Hamiltonian H(x, s, n x , n s ) as the starting point because 
our formalism requires a constraint which, with a given choice of variables, admits an 
intrinsic time, and at the level of this Hamiltonian there is no reason to prefer one definite 
sign for the non vanishing momentum n s . Anyway, it must be signaled that, because 
the reduced Hamiltonian rj/P is conserved, and in particular it does not vanish, in our 
interpretation there are no transitions from states on one sheet to states on the other sheet 
of the constraint surface, and therefore both points of view do not lead to an essential 
contradiction. 




2n s sign{n s )Q 




(27) 
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3 Canonical quantization 

3.1 Standard procedure 

In the literature we can find different solutions for the Taub universe. An interesting 
example within those which do not start from an explicit deparametrization is the solution 
found by Moncrief and Ryan || in the context of an analysis of the Bianchi type-IX 
universe with a rather general factor ordering of the Hamiltonian constraint [Q. In the 
case of the most trivial ordering they solved the Wheeler-DeWitt equation to obtain a 
wave function of the form 

P+) = jT duF{u)K lul (ie 20 "^) K 2w (~e 2 ^+) . (28) 

where K are the modified Bessel functions. In the particular case that F(u) = uj sinh(7ru;) 
they have shown that the wave function can be written in the form 

^(n,(3 + ) = R(n,f3 + )e~ s (29) 

with S a combination of exponentials in Q and j3 + . An important feature of this wave 
function is that for values of Q near the singularity (that is, the scale factor near zero) 
the probability is spread over all possible degrees of anisotropy given by j3 + , while for 
large values of the scale factor the probability is peaked around the isotropic Friedmann- 
Robert son- Walker closed universe. 

3.2 Schrodinger equation and boundary conditions 

A rather different approach beginning with the identification of a global phase time can be 
found in || , where the authors obtain a Schrodinger equation and its solutions are used to 
select a set of solutions of the Wheeler-DeWitt equation. They start from a Hamiltonian 
like that given in (f|), and solve the Wheeler-DeWitt equation 

/ d 2 d 2 1 4 \ 

(^-^-r 4l+ H^' !,)=0 (30) 

ll 



like Moncrief and Ryan to obtain a set of solutions of the form 

<Pu(,x,y) = 



a(uj)iu\e y ) + b{u)Ku\e y ) 



x 



c{u)I iui/2 {~e 2x ) +d(u)K iuj/2 (le 2x ' 
o o 



(31) 



with I and K the modified Bessel functions. Then they consider a change of variables 
analogous to ([]) but with only the minus sign, so that the momentum ir s is negative 
definite, and the time is then t — s; hence in (|7|) the first factor is positive definite, and 
the second one is a constraint linear in tt s = -n t which leads to the Schrodinger equation 



According to this interpretation, the contribution of the modified functions I iuJ /2 (|e 2x ) 
is discarded, because they diverge in the classicaly forbidden region associated to the 
exponential potential §e 4x ; the functions Jj a ,(|e 2/ ), instead, are not discarded, because in 
this picture the coordinate y is associated to the definition of time. 

As it was mentioned before, an interesting point to be noted is the fact that this time 
t coincides with the definition (|27|), including the ambiguity in its sign which is related 
with the existence of two sheets of the constraint surface. 

A point that should be remarked is the fact that the utilization of the Schrodinger 
equation in order to select the physical solutions of the Wheeler-De Witt equation has 
been possible because the potential in the square root does not depend on time. This 
ensures that the constraints fl3"0| ) and (|3"2|), which are equivalent at the classical level, 
are also equivalent in their quantum version as no additional terms arise associated with 
conmutators M. 
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3.3 Wheeler— DeWitt equation with extrinsic time 

Consequently with the above remarks we propose to solve a Wheeler-De Witt equation 
starting from a form of the Hamiltonian constraint such that a global phase time is easily 
identified among the canonical coordinates; hence, the resulting wave function has an 
evolutionary form and it could be understood as it is in ordinary quantum mechanics. 
The constraint (f|) allows to trivially define the time as 

t = —s sign(n s ). (33) 

We can make the usual substitution p k — ► — i-J^ to obtain the Wheeler-De Witt equa- 
tion 

' d 2 & 1 ^ 



\ dx 2 ds 2 9 
This equation admits the set of solutions 



e 4 * tf(x,s) = 0. (34) 



V u (x,s) = [a{uo)e lws + b{u))e~ iu 

x c{uo)I iu)/2 {\e 2x ) + d(uj)K iu}/2 (^e 2x ) 

j 

(35) 

where ±s is a global phase time. 

There are some points which deserve certain analysis: The first one is that we have 
chosen the real values for a>; the reason for this choice is the requirement to obtain the free 
particle solutions in the asymptotically free regime. The second point is that we have not 
discarded negative energy solutions; this would have been equivalent to select one sheet 
of the constraint surface. As we pointed in the context of path integral quantization, 
once we decide to work and give the results in terms of the variables (x, s, n x , n s ) there is 
no reason to choose one sign for the momentum tc s . Moreover, if we take the quadratic 
form of the constraint as an essential feature of gravitation, we should only admit a 
canonical transformation leading to a constraint equivalent to the original one; hence, 
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both signs of f\ must be considered simultaneously to ensure that the original constraint 
and that in terms of the new coordinates and momenta yield differential equations with 
the same number of solutions. One half of the solutions of our Wheeler-DeWitt equation 
correspond to those of the preceding section, which yielded from a Schrodinger equation. 
An important fact is that our procedure allows to obtain them without the necessity of 
defining a prescription for the square root operator, but only by choosing the trivial factor 
ordering. Finally, we should stress that although the original momenta cannot be avoided 
in the description, their role is restricted to appearing in the time variable, but not in the 
physical degree of freedom x. 

3.4 Another extrinsic time for the Taub universe 

On the other hand, in the case of the Taub universe a time in terms of the original variables 
can be found in a straightforward way by identifying the coordinate Q° conjugated to 
P = H. The Hamilton-Jacobi equation associated to the constraint H is 

The solution is clearly of the form 

W = W 1 (x) + W 2 (y) (37) 

where 

W x {x) = ±^J V3a 2 - e 4x dx 

with ± = sign(ir x ), and 

W 2 (y) = ± l -J v /3(a 2 -E)-4e 2 fdy 

with ± = signify). Matching the constants a and E to the new momenta P and P and 
following the procedure of Chapter 3 we have 



<3° 



" d 






(4/ 


_dP 





J Po=0 
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3P 2 - 4e 2 v 



3P 2 + V3P 2 -4e 2 ^ 



(3? 



with — for 7i y > and + for -n y < 0. Because on the constraint surface Pq = we have 



3tt 2 + (4/3)e 2 f = a 2 = P 2 then 7r y = ±i\/3P 2 - 4e 2 ^ and hence for both ir y > and 
< we obtain 

-o 1 / x /97r 2 + 4e 2 f + 37r y \ 

4y / 9vr 2 + 4e^ \ ^9vr 2 + 4e 2 f - 3^/ 

The gauge can be fixed by means of the canonical condition x = Q° ~ T{t) = with T 
a monotonic function. Thus, as 37r 2 + |e 2y = a 2 > 0, we can define an extrinsic time as 

t(n v ) = U\a\Q 
V3 

V3q; 2 — 3tt. 



i-^O 



In f vg + 3 M . (40) 



Now, if we go back to the original variables (Q, /3 + , 7Tq, 7r + ) we obtain 



t( 7 r„ 7 r + ) = lnf fg- ( ^ + 7r+) V (41) 

Again, it could be interesting to notice that the existence of different well defined global 
extrinsic times confirms that the coincidence of the time fl20l) and (|27|) shows the close 
relation existing between the path integral approach and the analysis of the reduction 
proposed in reference 0. 

4 Discussion 

The requirement of definition of a global phase time in terms of the momenta could sug- 
gest that we should abandon the idea of obtaining an amplitude for states characterized 
by the original coordinates. However, while a deparametrization in terms of the momenta 
may be completely valid at the classical level, it has been pointed by Barvinsky that at 
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the quantum level there is an obstacle which is peculiar of gravitation: There are basi- 
cally two representations for quantum operators, the coordinate representation and the 
representation in which the states are characterized by occupation numbers associated to 
given values of the momenta. The last one is appropriate when the theory under consid- 
eration allows for the existence of assimptotically free states associated to an adiabatic 
vanishing of interactions, so that a natural one-particle interpretation in terms of creation 
and annihilation operators exists. In quantum cosmology these assymptotic states do not, 
in general, exist. The suitable representation must be able to handle with essentially non 
linear and non polynomial interactions, and such a representation is a coordinate one. 
In this representation the quantum states are represented by wave functions in terms of 
the coordinates. The usual Dirac-Wheeler-DeWitt quantization with momentum oper- 
ators in the coordinate representation acting on *&(q) follows this line; but, as we have 
already observed, this formalism is devoided of a clear notion of time and evolution, un- 
less the potential is everywhere non null so that we can find a time among the canonical 
coordinates. 

Within this context, it would be interesting to explore the canonical transformations 
which lead to classical systems with a well defined intrinsic time as a systematic procedure 
of deparametrization. In fact, this it is one of the central aspects of the analysis presented 
in ||, where the canonical transformation (]?]) allows to obtain a system for which the 
Hamiltonian constraint admits the possibility to define an intrinsic time (i.e. in terms of 
the new variables). This scheme of work is justified by the observation above. 

Taking into account the importance of this proposal, here we have discussed the re- 
lation between the mechanism of reduction of the Taub cosmological model proposed 
by Catren and Ferraro and the path integral approach, finding that both lead to the 
identification of the same extrinsic time in terms of the original variables. We have also 
commented the way in which the ambiguity in the election of the constraint sheet mani- 
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festly appears in both frameworks. On the other hand, a secondary result of our work is 
the fact that this is the first application of the path integral approach presented in B to 
the Bianchi universes. 
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